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ABSTRACT 

We have simulated the formation of large-scale structure arising from COBE- 
normalized spectra computed by convolving a primordial double- inflation perturbation 
spectrum with the CDM transfer function. Due to the broken scale invariance ('BSI') 
characterizing the primordial perturbation spectrum, this model has less small-scale 
power than the (COBE-normalized) standard CDM model. The particle-mesh code 
(with 512 3 cells and 256 3 particles) includes a model for thermodynamic evolution of 
baryons in addition to the usual gravitational dynamics of dark matter. It provides 
an estimate of the local gas temperature. In particular, our galaxy-finding procedure 
seeks peaks in the distribution of gas that has cooled. It exploits the fact that "cold" 
particles trace visible matter better than average and thus provides a natural biasing 
mechanism. The basic picture of large-scale structure formation in the BSI model is 
the familiar hierarchical clustering scenario. We obtain particle in cell statistics, the 
galaxy correlation function, the cluster abundance and the cluster-cluster correlation 
function and statistics for large and small scale velocity fields. We also report here on 
a semi-quantitative study of the distribution of gas in different temperature ranges. 
Based on confrontation with observations and comparison with standard CDM, we 
conclude that the BSI scenario could represent a promising modification of the CDM 
picture capable of describing many details of large-scale structure formation. 

Key words: primordial power spectrum - CDM - large-scale structure - galaxy 
clustering - cosmological velocity fields 
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C<3 1 INTRODUCTION 1983, Mo et al. 1993). Other observational quantities include 

large-scale velocity fields in general, the line-of-sight distri- 
The search for a self-consistent model for the formation of bution of quasar absorption clouds, X-ray observations giv- 
large-scale structure capable of explaining the vast range in g information on the temperature and distribution of hot 
of available observational data poses one of the greatest gas, the observed network of filaments and "walls", the ex- 
challenges in theoretical astrophysics. In this paper, we will istence of voids and their observed properties, and evidence 
present strong evidence that a large body of observations is for the epoch of galaxy and cluster formation, 
consistent with a picture of structure formation based on a Prior to the analysis of the APM catalog and the COBE 
"double- inflationary" model (to be reviewed shortly). The measurement of the large-scale microwave anisotropy, the 
most directly relevant observations include the measured biased, as-normalized, flat CDM model (Blumenthal et al. 
large-scale anisotropy of the microwave background spec- 1982) was quite successful in predicting the observed hier- 
trum (Smoot etal. 1992, Gorski etal. 1994), power spec- archical network of filaments, pancakes and nodes and ob- 
tra analysis of galaxy catalogs (Park etal. 1992, Vogeley served galaxy clustering (White etal. 1987). Biasing, orig- 
et al. 1992, Fisher et al. 1993), count-in-cell analysis of galax- inally devised for explaining the enhanced correlation of 
ies (Efstathiou etal., 1990, Loveday etal. 1992), the galaxy Abell clusters (Kaiser 1984), also had the advantage of re- 
correlation function based on IRAS, CfA, and other cat- ducing small-scale velocity dispersions to reasonable levels, 
alogs (Loveday et al. 1992, Vogeley et al. 1992, Baugh & These advantages disappear with the COBE normalization 
Efstathiou 1993), the angular correlation function of APM of the CDM spectrum, because the antibiasing needed to ex- 
galaxies (Maddox et al. 1990) , abundance of clusters (Bah- plain the observed variances of counts in cells is quite difR- 
call & Cen, 1993)) as a function of mass, rms line-of-sight cult to reconcile with our observational knowledge of velocity 
peculiar relative velocity of galaxy pairs (Davis & Peebles, fields and with other tests. More precisely: The normaliza- 
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tion fixed by the COBE anisotropy measurement appears 
to imply excessive power in the small-scale (< 10fo -1 Mpc) 
regime of the primordial fluctuation spectrum P(k), which 
manifests itself in various conflicts with observations. Pre- 
dicted local velocity fluctuations seem to be quite in excess 
of measured relative pairwise line-of-sight projected veloc- 
ities (Davis & Peebles 1983), although the measurements 
may be sensitive to sampling effects (Mo et al. 1993) . 

COBE-normalized (flat) CDM simulations also predict 
much more pronounced and numerous bound structures 
such as massive galaxies and galaxy clusters than observed 
(Efstathiou et al. 1992). Relative to these massive structures, 
COBE-normalized CDM appears to imply fewer galaxies in 
filamentary structures than would be expected from obser- 
vations. In contrast, as we shall see, the double-inflationary 
models reproduce the observed filamentary structures rather 
easily and predict fewer massive galaxies and clusters. The 
different evolution of filamentary structures may be intu- 
itively understood within the Zel'dovich (1970) picture of 
anisotropic collapse of overdense regions leading to the for- 
mation of pancakes on a broad range of scales. In partic- 
ular, filaments are transitory structures: Due to the high 
normalisation, the COBE-normalized CDM is further de- 
veloped, and therefore it is plausible that cooled material 
in filaments will have more time to flow into large clumps 
(compare Doroshkevich et al. 1995) . 

On the other hand, "tilted" models (compare e.g. Cen 
et al. 1993) seem to have insufficient small-scale power: as 
a consequence, the predicted epoch of galaxy formation is 
difficult to reconcile with the observation of high-redshift 
objects. Moderate-scale (60ft~ 1 Mpc) power may also be too 
low, resulting in rather small streaming velocities. Mixed 
dark matter models (Davis et al. 1992, Klypin et al. 1993, 
Cen & Ostriker 1994) appear to agree better with a large 
body of observations, but again the late epoch of galaxy 
formation may be problematic (Mo & Miralda-Escude, 1994, 
Kauffmann & Chariot, 1994). 

Estimates of hot gas in clusters derived from X-ray 
observations seem to imply that a high percentage of the 
dynamical mass may be in the form of baryons. E.g., for 
the Coma cluster, White etal. (1994) compute that about 
30 % of the dynamical mass is in the form of hot gas. If 
the estimate for Qbh 2 w 0.013 obtained from nucleosynthe- 
sis is applicable in clusters (unless significant segregation 
of hot gas in clusters has occurred), then the total mat- 
ter density in the universe would be restricted to the range 
0.1 < f2 < 0.3. Indeed, low-mass CDM models (Kofman 
et al. 1993, Cen, Gnedin, & Ostriker 1993) with nonvanishing 
cosmological constant or low-mass and open CDM models 
(Kamionkowski et al. 1994) can explain the enhanced APM 
angular correlation function, probably the extended range 
of positive correlation of spatial galaxy correlation above 
30h~ Mpc, the cluster mass function, and the cluster-cluster 
correlation function; however, the epoch of galaxy and clus- 
ter formation appears to be rather early. In particular, clus- 
ters in open models tend to be highly developed and tightly 
bound, exhibiting less substructure than observed (Kaiser 
1991, Evrard etal., 1993). 

In addition to strictly observational constraints, the 
theoretical appeal of a model and its consistency with our 
knowledge of particle physics also play an important role in 
its acceptance. Standard CDM was very appealing because it 



involved only one true fit parameter (the amplitude of pri- 
mordial fluctuations) together with one phenomenological 
parameter (biasing). "Low-mass" models are theoretically 
less attractive because they require a nonvanishing cosmo- 
logical constant in order to achieve a flat universe, as pre- 
ferred in an inflationary scenario. Tilted models introduce 
one additional parameter; however, the simplest scenario as- 
suming an exponential inflaton potential is not theoretically 
well founded. Mixed dark matter models also involve only 
one additional parameter: (essentially) the neutrino mass 
(since their number density is fixed). The double-inflationary 
model requires two independent parameters. In Section 2, 
we review this model and describe fits to data in the lin- 
ear regime (Gottlober, Miiller & Starobinsky 1991, GMS91; 
Gottlober, Miicket & Starobinsky, 1994, GMS94) and in the 
Press-Schechter theory (Miiller 1994a) resulting in the par- 
ticular "BSI" (Broken Scale Invariance) spectrum studied in 
this paper. 

Encouraged by the success of the BSI model in the lin- 
ear regime, we now present the results of numerical n-body 
simulations in order to test observational results requiring 
accurate predictions in the nonlinear regime. Now, one of 
the subtleties of the nonlinear regime is that our incom- 
plete understanding of physical processes at relevant length 
scales causes theoretical uncertainties that may be compara- 
ble to observational uncertainties (cp. for example, Ostriker 
1993). Complicated, nonlinear feedback mechanisms govern 
the dynamics at scales below our numerical resolution. In 
particular, an adequate theory of galaxy formation (and per- 
haps evolution) together with significantly improved com- 
putational resources would obviate the necessity for intro- 
ducing simple phenomenological concepts such as biasing. 
However, an accurate description of galaxy formation re- 
quires a robust treatment of gas dynamics over vast ranges 
of scales, densities, and other conditions; dynamic computa- 
tion of heating and cooling, including thermal instability and 
ionization; star formation, supernovae; proper handling of 
the chemical evolution of a multi-component medium — to 
name a few considerations. Significant progress in applying 
a hydrodynamic approach to cosmology has been achieved 
by Katz & Gunn (1991), Navarro & Benz (1991), Cen & Os- 
triker (1992), Katz, Hernquist & Weinberg (1992), and by 
Steinmetz & Miiller (1995). In particular, a hydrodynami- 
cal description is essential for studying the formation and 
internal dynamics of clusters (Evrard 1990). A compara- 
tive review containing additional references is now available 
(Kang et al. 1994). An alternative hydrodynamical approach 
emphasizing thermal instability and the importance of feed- 
back mechanisms such as supernovae has been studied by 
Klypin, Kates & Khokhlov (1992; KKK92). Hydrodynamic 
simulations are clearly superior for obtaining reliable small- 
scale predictions from models of large-scale structure forma- 
tion. However, a lack of computer resources still precludes 
their general use for high-resolution cosmological applica- 
tions. 

The numerical particle-mesh code used here (Kates 
etal. 1991, KKK91; Klypin & Kates 1992, KK) includes a 
model for thermodynamic evolution of baryons in addition 
to the usual gravitational dynamics of dark matter (see Sec- 
tion 3). The code provides an estimate of the (suitably aver- 
aged) local gas temperature without the enormous computa- 
tional cost of a full hydrodynamical description. A resolution 
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of 512 3 cells (256 3 particles) was achieved. High resolution 
is crucial in the present study in order to include the effects 
of both enhanced large-scale power and modest (compared 
to CDM) small-scale power in the same simulation. 

As discussed in Section ^, considerable preliminary test- 
ing using both linear analysis (compare GMS94) and moder- 
ate (256 3 cells) resolution (compare Gottlober 1994, Miiller 
1994b, Amendola et al. 1995) was first carried out to obtain 

ftimal estimates for the parameters kb r and A (see Section 
). These preliminary investigations also provided rough 
estimates of the sensitivity of our results to "fine-tuning" 
of parameters (see Conclusions), box size, etc. The power 
spectrum used in the highest-resolution simulations and its 
implementation in the numerical simulations will be given 
in Section |ij 

Simulation of the coupling of thermodynamic evolution 
to gravitational dynamics provides us with important indi- 
rect information concerning observable quantities; see Sec- 
tion ^j. In Section |f| we present our results on the evolu- 
tion of the mass distribution and, in Section |(| its relation- 
ship to thermodynamic evolution. Of particular interest are 
the particle in cell fluctuations, the evolution of clustering 
with redshift and its dependence on the spectrum, redshift- 
dependence of thermodynamics and percentage of cooled 
particles, the relation between density peaks and hot and 
cold particles. 

In Sections (?] and ^ we discuss the results of galaxy iden- 
tification and galaxy clustering: galaxy mass function versus 
Press-Schechter theory; dependence on spectrum and red- 
shift; galaxy correlation function versus correlation function 
of density peaks; clustering and other properties of galaxy 
clusters. The characteristic features of large- and small-scale 
velocity fields for BSI spectra are presented in Section ^. In 
Section lot we draw our conclusions. 



2 BSI POWER SPECTRUM 

Standard inflation produces primordial adiabatic density 
perturbations P(k) oc k n with a spectral index n = 1 plus 
small logarithmic corrections which stem from the slowly 
changing horizon length during inflation, i.e. they depend 
on the inflaton potential. We propose an early cosmological 
evolution which produces a power spectrum with a break 
at a characteristic scale by introducing more than one effec- 
tive field responsible for inflation (Starobinsky, 1985). We 
use the combined action of renormalization corrections and 
a massive scalar field as the source of a non-flat primordial 
perturbation spectrum (GMS91). This model decouples the 
clustering properties at large and small scales. It requires 
two new parameters in addition to the amplitude: First, the 
ratio of the masses A ~ m/(6.5A/) appearing in the La- 
grangian 

this ratio essentially determines the ratio of power at small 
(galactic) scales to that at very large cosmic scales. Second, 
the epoch of the transition between inflationary phases is 
related to the initial energy density m 2 cf>Q of the scalar field. 
This epoch determines a characteristic length scale ^ r = 
2ir/k-fo I where the shape of the spectrum changes. 




0.001 0.010 0.100 1.000 10.000 
k [h Mpc~'] 



Figure 1. BSI power spectrum (solid line) compared with stan- 
dard CDM (dashed line), a ACDM model (dotted line), and a 
MDM spectrum (dash-dotted line). 

In GMS94, linear theory predictions of double- 
inflationary models were compared with observational con- 
straints. Assuming COBE normalization, the tests applied 
were as follows: "counts in cells" of the IRAS and APM sur- 
veys, the APM galaxy angular correlation function, bulk- 
flow peculiar velocities, the Mach number test, and quasar 
abundance. The BSI models were shown to be in good agree- 
ment with all of these tests for the parameter regime defined 
by 2 < A < 4 and 0.5/i _1 Mpc< k7^ < 5/i _1 Mpc, the lat- 
ter corresponding to length scales iu r ~ (3 — 30)ft _1 Mpc. 
Although some freedom is permitted by these constraints in 
the choice of k^ r (or ^ r ), the best fit seems to be given by 
fcr 1 = 1.5ft" 1 Mpc, A = 3. These values also led to reason- 
able effective biasing parameters in preliminary test simula- 
tions. In all simulation results reported here, the terminol- 
ogy "BSI spectrum" refers to these values of the parameters 
unless otherwise stated. 

After convolution with the standard CDM transfer 
function of Bond & Efstathiou (1984), taking Q = 1 and 
h = 0.5, the power spectrum is used to generate the initial 
condition of the numerical simulations described in the next 
section. Fig. |l| plots the transformed BSI spectrum in com- 
parison with the power spectra of the standard CDM model, 
a ACDM (A = A/3H 2 = 0.8, = 0.2) and an MDM model 

( n CDM = °- 7 > Q " = °- 2 ' Qb = al )- The standard CDM 
spectrum has much more power at all the relevant scales 
for galaxy formation, while BSI mostly resembles MDM at 
large scales, but has more power on small scales, almost as 
high as the ACDM model. An analytical fit to the primordial 
potential spectrum is given by 

, 3(6 , H _ / 4.2 x 10- 6 [log {k s /k)]°- e + 4.7 x lO" 6 (k < k s ) 
[ ' ~\ 9.4 x 10- 8 log (k f /k) (k > k s ) 

with k B = (27r/24)/iMpc"\ k f = e 56 /iMpc _1 . However, for 
a proper representation of the transition regime, the numer- 
ically tabulated spectrum is preferable. 

In Fig. ^| we show the predicted multipole moments of 
the anisotropies of CMB fluctuations calculated following 
Gottlober & Miicket (1991) in comparison with new mea- 
surements on the angular scale 2 < I < 1200. Our spec- 



© 1995 RAS, MNRAS 000, §-?? 



4 R. Kates, V. Miiller, S. Gottldber, J. P. Miicket, J. 



Retzlaff 





0.01 



BSI-25 



0.10 1.00 10.00 100.00 

k [h Mpc~' ] 



Figure 2. Comparison of the cosmic background fluctuation mul- 
tipolcs of BSI (solid line) and standard CDM (dashed line) with 
a series of experiments. 

tra are normalized with the 10°-variance of the fluctua- 
tions g t = (30 ± 7.5) /iK/2.735 K from the COBE first-year 
data (Smoot et al. 1992). The triangle denotes the first-year 
COBE result Q rms _ps = (16.7 ± 4)/xK (for the spectral 
index n = 1) which is approximately equivalent to our nor- 
malization. The other experimental data (courtesy of B. Ra- 
tra, from Bond 1994) are from left to right COBE, FIRS, 
Tenerife, SP91, SK93, the lower end of the PYTHON error, 
ARGO, MAX1 and MAX2, full and source free MSAM2 
and MSAM3, and the upper limits of WD and OVRO (for 
details see Bond 1994). The reanalysis of the COBE second- 
year data by Gorski etal. (1994) implies an about 25 per 
cent higher normalization of the spectra (Gottlober 1994). 
Consequences of such an increase will be discussed in the 
conclusions. 



3 NUMERICAL REALIZATION 

Initial fluctuations corresponding to BSI and standard CDM 
power spectra were generated as realizations of a Gaussian 
random field. Positions and velocities were assigned accord- 
ing to the Zel'dovich approximation at redshift z — 25 (al- 
though for larger boxes the Zel'dovich approximation could 
have been continued until a later epoch). The Euler-Poisson 
system describing the evolution of self-gravitating, collision- 
less matter was then evolved from z = 25 to z — using the 
particle-mesh (PM) code described in KKK91 and KK. 

To cover a large range of the spectra, we performed 
a series of simulations with different box sizes (see Table 
1). Asterisks denote high resolution simulations(512 3 cells 
and 256 3 particles); the remainder of the simulations were 
performed using 256 3 cells and 128 3 particles. The two 
high-resolution BSI simulations were performed in boxes of 
200/i _1 Mpc and 25/i _1 Mpc. BSI-200* is intended for anal- 
ysis of the large-scale matter distribution, while BSI-25* 
should provide good working accuracy on galactic scales. 
In order to highlight differences between BSI and standard 
CDM, a series of simulations with half of this resolution 
were performed using the same random seed in boxes of 
size 25/!~ 1 Mpc, 75/i _1 Mpc, 200/i" 1 Mpc, and 500h -1 Mpc. 



Figure 3. Range of spectra realized in the simulations BSI-25 - 
500 (horizontal lines indicate fc mm to fcmax); also shown is the 
linear theory BSI spectrum (dash-dotted line) and a reconstruc- 
tion assembled from the BSI simulations at z = (full line). 

'grid Si ves the size of the particles used in the cloud-in-cell 
mass assignment scheme as well as the formal resolution of 
the grid on which the gravitational force is calculated. The 
Nyquist wavelength 27r/fcmax, which depends on the number 
of particles used in the simulation, is four times larger than 
'grid' ^ corresponds to an upper limit for the realization of 
initial perturbations in k-space. Mb+ d m is the total mass of 
a particle. Finally, we compare the grid variances of density 
contrast as and peculiar velocities a v of the initial realiza- 
tion of the Gaussian random field with the expectations of 
the linear power spectrum, 

fc max ^ fc max 

*« = 2^ / dfcfc2p ( fc )> ^=<|l / dkP(k), (3) 
fe min fc min 

where the numbers are given at z = transformed according 
to the linear theory growth law. They describe the quality 
of the realizations of the power spectra in the different sim- 
ulation boxes. 

On the scale of 8 /i _1 Mpc, one often refers to the mea- 
sured unit variance of galaxy counts in spheres as deter- 
mined from the CfA-catalog ((Davis & Peebles, 1983). Using 
the first integral of Eq. 3 with a top-hat window of radius 
8 /i _1 Mpc, we infer as = 0.46 for the BSI power spectrum. 
Therefore we predict a (linear) bias of galaxies with respect 
to dark matter of b ~ 2.2. For the standard CDM model 
we have on the other hand erf = 1.12, i.e. there the galaxies 
should be slightly antibiased. 



4 THERMODYNAMIC ESTIMATES OF GAS 
TEMPERATURES 

As described in KKK91 and KK, estimates of the average 
values of the local dark matter density and velocity and of 
the local baryon temperature can be obtained from simula- 
tions of large-scale structure (without direct simulation of 
hydrodynamics) by following the thermal history of the gas 
while imagining that, smeared out over a sufficiently coarse 
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Table 1. Parameters of the simulations 



Simulation 'grid %+DM a S &5 u v &v 

[fo-iMpc] [M Q ] [km s" 1 ] [km s" 1 ] 



BSI-200* 


0.39 


2.6 


X 


10 11 


2.28 


2.25 


490 


483 


BSI-25* 


0.05 


5.2 


X 


10 08 


4.60 


4.69 


223 


219 


BSI-500 


1.95 


3.3 


X 


10 13 


0.94 


0.93 


512 


508 


CDM-500 


1.95 


3.3 


X 


10 13 


3.02 


2.99 


1049 


1040 


BSI-200 


0.78 


2.1 


X 


10 12 


1.62 


1.60 


516 


479 


CDM-200 


0.78 


2.1 


X 


10 12 


5.43 


5.36 


1103 


1068 


BSI-75 


0.29 


1.1 


X 


10 11 


2.59 


2.54 


441 


376 


CDM-75 


0.29 


1.1 


X 


10 11 


8.79 


8.65 


1064 


979 


BSI-25 


0.10 


4.1 


X 


10 09 


3.85 


3.76 


269 


219 


CDM-25 


0.10 


4.1 


X 


10 09 


13.1 


12.8 


857 


735 



scale, baryons are transported with the dark matter; i.e., 
pgas = P^b' w here flu = 0.1 is the background fraction of 
baryons. The practical advantage of this picture is that the 
thermal history of the gas reduces to integration of one or 
more ordinary differential equations along known trajecto- 
ries (those of the dark matter). 

This description is certainly unproblematic before the 
formation of the first shocks, when the medium is still cold 
and fluctuations are small (we put Tgas = at 2 s tart = 25). 
As perturbations grow, eventually the first objects start to 
collapse, producing caustics in the dark matter and shocks in 
the gas. As simple pancake models show (Shapiro & Struck- 
Marcell 1985), shocks occur close to caustics. 

The label "shocked" may be assigned to a particle in one 
of two ways: First, particles for which the Jacobian deter- 
minant of the transformation from Lagrangian to Eulerian 
coordinates is negative are classified as shocked. (Particle 
labels (i,j,k) are the numerical realization of Lagrangian 
coordinates. The position assigned to a particle is the nu- 
merical representation of its Eulerian coordinates.) The Ja- 
cobian is determined for each particle by considering the 
volume element constructed from the position vectors to its 
nearest Lagrangian neighbors. Second, a particle inherits the 
label "shocked" by passing through a shocked region. More 
precisely, if particle (i,j,k) is in cell (I,J,K), then the la- 
bel shocked is assigned to particle (i,j,k) if the density of 
previously shocked gas is nonzero in cell (7, J, K) and all 
neighboring cells. 

At a shock, the temperature acquired by gas particles 
with velocity v is given by 

kT « ti M m H {v~U) 2 /3, (4) 

where U is the local velocity determined by interpolation of 
the velocity field onto a coarse grid of twice the usual cell 
size. (If it should happen that too few particles are present 
to determine U, the temperature assignment is simply post- 
poned.) This estimate is relatively insensitive to small er- 
rors in the position of the shock. Here, mn is the mass 
of hydrogen, and /j,m is the molecular weight per particle, 
Mm ~ (n H + 4n H e)/(2n H + 3n He ) ~ 0.6. 

When a particle crosses a shock and is assigned a tem- 
perature, we start to integrate the energy equation along the 
trajectory of the particle: 



where /j,h is the molecular weight per hydrogen atom, fiH ~ 
(n H + 4n H e)/(nH) ~ 1-4, pgas = pflu, = n H ^Hm H , and 
7 = 5/3. Here, A ra( j represents radiative losses in the hot 
plasma with assumed primordial abundances, and AQ om p 
is the cooling rate due to Compton scattering. To estimate 
A ra( j, we used analytic fits as in KK and KKK91 to the 
cooling curves given in Fall & Rees (1985). The medium is 
treated as optically thin and in collisional equilibrium. 

Computation of the increment in T due to the (La- 
grangian) time derivative of nu in Eq. |^ requires a knowl- 
edge of uh at the present and at the previous timestep. For 
each timestep, we construct the "gas density" on the grid, 
defined by counting only contributions of shocked (hot) par- 
ticles. (For the dark matter density, one of course counts all 
particles). For the evolution of Eq. || the present value of 
nu for Eq. [B| is computed simply by interpolation to the po- 
sition of the particle and then stored for use at the following 
timestep. 

Now, when material cools below 10 4 K, stars will form, 
producing luminous matter. By assigning the label "cooled" 
to the particles with T < 10 4 K, we have a measure for the 
amount of visible matter. However, we know that the effi- 
ciency for the conversion of gas to stars is low. One nonlin- 
ear feedback mechanism limiting star formation could be the 
energy provided to the medium by supernova explosions. A 
second (also nonlinear) mechanism inhibiting star formation 
is ionization and heating due to ambient ultraviolet radia- 
tion. As emphasized by Efstathiou (1992), a similar mecha- 
nism could be responsible for suppressing the formation of 
dwarf galaxies. It should be emphasized that some mecha- 
nism preventing conversion of the gas to luminous matter is 
necessary in BSI models as well as in CDM and in the pre- 
COBE version of CDM: otherwise, all the gas would simply 
collapse, forming stars and ending up in globular clusters 
long before the formation of large objects. 

In the present realization of the code, we take these 
effects into account in a crude and purely local manner: 
we "reheat" those particles to the temperature T renea t — 
5 x 10 4 K which have cooled to below 10 4 K with a proba- 
bility -Preheat and assign to the remaining particles the la- 
bel "cooled." Some of the reheated particles later will attain 
high temperatures, provided they enter collapsing regions of 
high density. Reheated particles can simply cool again and 
turn into "visible" matter. However, once cooled, a particle 
cannot be reshocked. The value f ren eat was held constant 
at 85 per cent, the value estimated in KK for a 50 /t -1 Mpc 
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grid and a as normalized CDM spectrum to give a reason- 
able fraction of hot gas in clumps of galactic mass (about 
10 per cent). The fraction of cooled gas (88 per cent) for a 
cell size of about 0.05/i -1 Mpc in the highest resolution BSI 
simulation seems to be reasonable. The distribution of par- 
ticles of different temperature ranges, variation with cell size 
and spectrum, and related considerations will be discussed 
below in Section ^. 

We recall that the mean density of stars (and other 
components of the baryonic matter) in some comoving cell 
depends in principle on the entire past history of the cell and 
not just on its density at some particular epoch. As discussed 
in KKK91 and KK, comparison of the density distribution 
of "cold" material - which is related to the density of visible 
matter - with the dark-matter density thus yields informa- 
tion useful for testing cooling as a physical mechanism for 
"biased galaxy formation." 

A reasonable upper bound for the aforementioned 
"coarse graining" (i.e., error in predicting positions of 
baryons) is the local sound velocity integrated over a Hubble 
time (e.g., for T « 5 x 10 6 K it is 2 or 3fr _1 Mpc). However, in 
tests - leaving aside conditions prevailing in clusters - for a 
box of 75 ft _1 Mpc (cell length 200 /i _1 kpc), the error seems 
to be of the order of the cell resolution only. The explanation 
seems to be that both components move in the same poten- 
tial well. The dark matter spends a large fraction of its time 
at about the same radius as the gas, because they originally 
had the same kinetic energy. Due to mixing, the dynam- 
ics of dark matter particles bears some resemblance to the 
dynamics of a gas whose temperature corresponds to the lo- 
cal velocity dispersion. Incidentally, from experiments with 
truncated spectra in KKK91 and from theoretical consider- 
ations, we know that temperature predictions are somewhat 
sensitive to the portion of the spectrum actually simulated. 
This fact implies that the same spectrum simulated in boxes 
of varying sizes can yield varying estimates of temperature 
distributions. 

Results of a ID-pancake test presented in KKK91 show 
good agreement of our model with hydrodynamical simula- 
tions (20 per cent of gas cooled, position of cooling front and 
shock wave). Proper inclusion of hydrodynamics and the ef- 
fects of thermal instability would be expected to make a sig- 
nificant difference in at least two situations: First, when gas 
starts to cool efficiently - this happens inside dense regions 
and/or if the temperature becomes too low (T < 2 x 10 5 K). 
Second, when gas undergoes secondary shocking - for ex- 
ample, in collapse to objects with masses smaller than a 
galactic mass (for more detailed comparison, see, KKK92). 
Thus, the code does not properly treat the internal regions 
of galaxies and clusters, but it would be expected to give a 
reasonable approximation for the temperature and distribu- 
tion of gas which leaves voids and is trapped in the potential 
wells of superclusters and filaments. Comparison with Cen & 
Ostriker (1992) supports this expectation: the baryonic and 
dark matter distributions look remarkably similar; most of 
the gas in regions with density p > 10(p) has a tempera- 
ture between 10 6 K and 10 7 K, similarly to what was found 
in KK. Comparison with hydrodynamic simulations of peri- 
odic disturbances reported in KKK92 also lead to reasonable 
qualitative agreement in this regime. 
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Figure 4. Comparison of the variances of counts in cells of galax- 
ies of the Stromlo-APM (full circles) and IRAS (triangles) survey 
with the variances of the dark matter particles in the simulations. 
The BSI curve (solid line) composed from the simulations with 
different box sizes fits the data quite well. The dashed line corre- 
sponds to the CDM simulations. Dashed-dotted and dotted lines 
are the linear theory predictions of BSI and CDM, respectively. 



5 COUNT-IN-CELL ANALYSIS FOR DARK 
MATTER 

One useful method of characterizing the dynamics of clus- 
tering is to study the statistics of counts in cells. It is a 
robust measure for distinguishing different power spectra. 
For example, the rms fluctuations of simulated dark matter 
particles in a sphere of radius 8/i _1 Mpc are often used to 
determine the linear biasing factor. The fluctuations on the 
scale of a galactic halo provide an estimate of the threshold 
density contrast required in galaxy identification algorithms. 
Comparison of the simulated count fluctuations as a function 
of radius with rms fluctuations observed in the IRAS cata- 
log (Efstathiou et al. 1990a) and the Stromlo-APM survey 
(Loveday et al. 1992) provides an important test of models. 

In Fig. ^|the curves representing COBE-normalized BSI 
have been shifted vertically by a "biasing" factor 6ggj = 2 in 
order to normalize fluctuations to observations at 8fa _1 Mpc, 
while the CDM curves had b^DM = °- 9 - CDM fluctua- 
tions are too low (significantly outside the 1-cr error range 
on the low side) on scales (15 — 50)/i -1 Mpc, after which the 
measurement errors increase. In contrast, the interpolation 
curves in the range of BSI simulations fit the data quite 
well. The good fit beyond about 15/i _1 Mpc is consistent 
with the predictions of linear analysis (GMS94). The slight 
differences at small scales are due to the difference between 
linear and nonlinear evolution. The variances of cell counts 
are calculated in real space. At the scales studied the linear 
theory predicts a constant amplification factor which can be 
incorporated in the value of the bias parameter. The mea- 
sured variances provide a good test of the slope of the BSI 
power spectrum. In contrast, the CDM counts in cells are 
inconsistent with observations at about the 2-a level. 

The redshift dependence of the probability distribution 
of particle-in-cell counts delivers an estimate of the non- 
linearity of clustering. Using the counts on a 128 3 grid of 
the 200 fo -1 Mpc-simulations, we compare BSI and CDM 
simulations in Fig. B. Over four decades in the probability 
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Figure 5. Particle-in-cell distribution of BSI-200 (full lines) 
and CDM-200 (dashed lines) simulations. The curves (z = 
2,1.5,1,0.5,0) show the dependence of the probability distribu- 
tion of grid cells on the number of particles in the cell; the curves 
are steeper at higher redshift. For reference, a curve of slope —3 is 
drawn, corresponding to an evolved stage of pancake formation. 



it is remarkably well described by a power law (f>(n) oc n~ 3 , 
which is a prediction of the pancake model (Kofman etal. 
1993). The CDM models at z — are less steep, correspond- 
ing to more big clumps in the matter distribution. The CDM 
models at about z = 2 had a similar slope to the BSI model 
at z = 0. The subsequent evolution to a less rapid fall off is a 
quantitative measure of the substantially increased cluster- 
ing of CDM compared to BSI models. The large dynamical 
range of the dark-matter probability distribution as com- 
puted using the simulation cells is more of theoretical inter- 
est, and it cannot be compared directly with the one-point 
distribution function estimated from galaxy surveys. 



6 GAS TEMPERATURES: STATISTICS AND 
SPATIAL STRUCTURE 

As explained in Section ^, the "temperature" assigned to 
each particle depends on its complete history, including the 
epoch and conditions under which it was shocked and the 
environment (in particular the local density) at all interme- 
diate times since shocking. At a particular time, say z = 0, 
the gas temperature distribution thus contains more infor- 
mation than the density and velocity distributions alone. 

In interpreting the numerically computed gas tempera- 
ture distribution, one should keep in mind the idealizations 
and limitations involved. As discussed above, a proper hy- 
drodynamic treatment of the gas would be expected to make 
a systematic difference in regions with secondary shocking 
and/or very high density, in particular inside clusters and 
galaxies. However, any treatment of the gas, even including 
hydrodynamics, will contain numerical inaccuracies (result- 
ing both from limited resolution and from systematic errors) 
as well as theoretical uncertainties due to incomplete mod- 
elling. Bearing all of these errors and uncertainties in mind, 
we take a brief look at what can be learned from particle 
temperature statistics. 

In Fig. ^, the percentage of particles in selected temper- 
ature ranges at z = for BSI and CDM simulations is plot- 
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Figure 6. Percentage of particles in different temperature ranges 
as a function of cell size at redshift 2 = 0. Plus signs indicate 
CDM simulations, asterisks BSI, and diamonds high-resolution 
(512 3 cells) BSI simulations, a) never shocked, b) cooled, c) 
temperature range 10 4 < T < 5 X 10 4 , d) temperature range 
10 5 < T < 1.5 X 10 6 , e) temperature range T > 10 7 . 



ted against cell size, which for all simulations shown is in- 
versely proportional to the highest resolved wavenumber. In- 
deed, studies of truncated spectra in high-resolution simula- 
tions (KKK91) indicate that the wavenumber limit is crucial 
to temperature statistics, i.e., qualitatively similar trends 
would be observed even if the spatial resolution were to be 
improved while keeping the wavenumber limit constant. (In- 
cidentally, percentages do not add up to one, because the 
temperature ranges shown are not exhaustive.) Plus signs 
indicate CDM simulations, asterisks and diamonds BSI; di- 
amonds indicate the high-resolution simulations, in which 
the ratio of highest to lowest resolved wavenumber is twice 
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Figure 7. Comparison of slices through 200/i 1 Mpc BSI (a) and 
CDM (b) simulations in the temperature range 10 4 K < T < 
5 X 10 4 K. 



as large. One expects the (limited) spectral range to have 
an influence on gas temperature statistics, and this expec- 
tation is borne out in Fig. |(| The trend with respect to cell 
size is different in different temperature ranges. However, the 
most striking regularity is that in all temperature regimes, 
a smooth (in most cases monotonic) trend is apparent upon 
plotting percentage against cell size (equivalently here: high- 
est resolved wavenumber) rather than 60a; size (equivalently 
here: lowest resolved wavenumber). A plot with respect to 
box size would be obtained by shifting the high-resolution 
points (diamonds) by one point to the right, which in all 
cases would move them off an otherwise smooth curve. 

Fig. [?] illustrates a typical slice in the range 10 4 K 
< T < 5 x 10 4 K for both BSI and CDM simulations. (Note 
that the random seed used in BSI and CDM simulations was 
the same, resulting in the same phase relations for each per- 
turbation mode.) Slices are about 12/i _1 Mpc thick. About 
9000 particles are plotted in each slice (all of the particles 
in the CDM slice and about 40 per cent of the particles in 
the BSI slice). In Fig. [s| spatial distributions for BSI simu- 
lations are illustrated in the same thin slice for particles in 
the ranges 10 5 K < T < 1.5 x 10 6 K, T > 10 7 K, cooled 
(T < 10 4 K), and never shocked. The different temperature 
ranges highlight quite different features of the distribution 
in a way that would not be possible using simple peak statis- 
tics. 

We now consider various temperature ranges in detail: 
Never shocked: This percentage is a measure of the frac- 
tion of particles not in pancakes. The fraction increases 
monotonically with cell size for both BSI and CDM (Fig. 
^). The trend may be attributable to the fact that for larger 
cell size, the portion of the spectrum actually simulated is 
shifted toward smaller wavenumbers, with correspondingly 
smaller values of k 3 P(k), resulting in longer timescales for 
formation of structures such as pancakes. The number of un- 
shocked particles is consistently lower by a factor of about 
three for CDM than for BSI, directly reflecting the signif- 
icantly enhanced small and medium-scale power of CDM 
compared to BSI. In Fig. ^a, the spatial distribution of un- 
shocked particles is shown for a 25/i -1 Mpc high-resolution 
simulation. The distribution shows practically no structure, 
aside from being absent in regions of high density, where 
all of the particles are shocked. Visual inspection of slices 
confirms one's expectation that material in voids consists 
mainly of unshocked particles. 



Cooled: From plausibility arguments one expects the distri- 
bution of "cooled" particles to trace the galaxy (halo) distri- 
bution more accurately than does the unbiased dark-matter 
distribution. (A higher probability of cold gas obviously fa- 
vors star formation.) The results of KKK91 in the context 
of high-resolution, 2D simulations lend additional support 
to this expectation: There, halos were found with high con- 
fidence by a ("friends-to- friends") clustering algorithm ap- 
plied to all particles. The percentage of cold particles inside 
halos was significantly higher than outside. Moreover, the 
epoch of galaxy formation (estimated by cluster analysis at 
various redshifts) agreed with the average redshift Zg a j at 
which the cold particles in galaxies "cooled." In three dimen- 
sions, we exploit the affinity of cooled particles for galaxies 
in our galaxy-finding routines: As discussed in Section [j], the 
routines begin by searching for peaks in the cooled-particle 
distribution. 

As seen in Fig. ^b, the percentage of cooled particles 
at z — varies strongly with cell length (more precisely, 
with Nyquist frequency), mainly because power at smaller 
length scales leads to early formation of small pancakes and, 
in spite of "reheating," subsequent cooling. (The redshift 
at cooling is recorded for each cooled particle.) Consider a 
structure of characteristic size L. Smaller structures are as- 
sociated on the average with smaller characteristic initial 
temperatures (oc L 2 , since typical relative velocities in Eq. ^ 
scale with HqL). The CDM simulations have more small- 
scale power and therefore consistently fewer unshocked par- 
ticles and more cooled particles. In KKK91, it was shown 
that truncation of the spectrum below the Nyquist frequency 
leads to significant reduction of the number of cooled par- 
ticles (keeping spatial resolution constant). Thus, as in the 
case of unshocked particles, the variations with "cell size" 
seen here can probably be attributed to changes in the range 
of the spectrum actually simulated. Filaments and voids are 
evident in the cold particle distribution; the hierarchy of fil- 
ament separations and void sizes give the visual impression 
of a broad distribution of characteristic length scales. 
Particles in the temperature range (1 — 5) x 10 4 K: 
This temperature range would tend to be associated with 
Lyman-a clouds. However, a more realistic treatment of the 
effect of the ultraviolet background on the matter (including 
ionization and heating) would be required to draw quantita- 
tive conclusions concerning Lyman-a clouds. Such a treat- 
ment is discussed by (Miicket etal. 1995). Nevertheless, we 
may obtain some hints as to what to expect from a more 
realistic treatment by examining the trend in Figs. ^p. The 
fraction grows monotonically both in BSI and in CDM with 
cell scale, except for the 500ft _1 MpcBSI simulation, with 
CDM consistently lower. A reasonable interpretation may 
be that earlier pancake formation and efficient cooling at 
moderate redshifts in the present code deplete the reservoir 
of particles (especially those near galaxies) that would other- 
wise be candidates for this temperature range. The turnover 
of the BSI curve at large scales is not surprising, since at 
the largest scales, the reservoir of shocked particles sets the 
limit. (In any case, quantitative studies of Lyman-a clouds 
would require resolution of about 50 — 100 kpc.) 

With proper treatment of ionization and heating, the 
relative fraction of cold particles would be expected to drop 
significantly in favor of particles in the range of tempera- 
tures (1 — 5) x 10 4 K, except in very dense regions. Upon 
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examination of numerous slices, there is the visual impres- 
sion that particles in the temperature range (1 — 5) x 10 4 K 
at least roughly trace the distribution of cold particles at 
2 = 0. This may be a hint that a population of Lyman-Q 
clouds could be associated with galaxies (cp. Petitjean etal. 
1995). 

Comparing BSI and CDM slices for the temperature 
range considered in Fig. 7, there is a strong impression of 
thinner, better defined filaments and emptier voids in CDM 
than in BSI simulations. (Note that the total number of 
particles is the same.) This impression is consistent with 
the interpretation that the standard CDM looks like the 
BSI model would look like if it were evolved to a larger 
fluctuation amplitude. 

Particles in the temperature range 10 5 K < T < 

1.5 x 10 6 K: As explained in KK, there are strong theoretical 
and observational reasons for an association between gas in 
this "warm" temperature range and filamentary structures. 
In Fig. ^c, it is evident that warm particles trace filaments 
quite closely in BSI. The same qualitative picture is found 
in all boxes. Comparison of CDM and BSI simulations in 
this warm regime shows little qualitative difference in fila- 
mentary structures, except that CDM consistently leads to 
somewhat more gas in the warm regime than BSI. This oc- 
curs despite CDM's enhanced accumulation of material in 
clusters, which would tend to transfer gas from the warm to 
the hot regime (see discussion of hot gas below). The trend 
toward more warm gas at larger simulations may be related 
both to the increased reservoir of shocked, but not yet cooled 
particles (an essentially small-scale effect) and perhaps the 
influence of large-scale power in the portion of the spectrum 
actually simulated. This second cause is consistent with a 
rise just visible in the BSI curve at cell-size 0.39/i _1 Mpc 
(high-resolution simulation) . 

Particles in the temperature range T > 10 7 K: The 

"hot" gas represented by these particles is generally associ- 
ated with the deep potential wells of clusters, and therefore 
it is no surprise that in Fig. 8d hot gas traces clusters (also 
verified by comparison with galaxy catalogs). The percent- 
age of "hot" particles with T > 10 7 K, which are nearly 
all associated with clusters or large galaxy groups, is consis- 
tently an order of magnitude higher for CDM than for BSI, 
the trend becoming even stronger at larger cell size. This 
trend is consistent with the general conclusion that CDM 
causes more and larger massive structures (also evident in 
slices). Moreover, these structures seem to be farther evolved 
in CDM, so that deep potential wells would also tend to be 
broader. Broader cluster potential wells would explain why 
the disparity between CDM and BSI grows at larger cell 
lengths (Fig. |^e) . Cluster potential wells are still fairly well 
resolved at the poorest resolution in CDM, but not in BSI. 



7 GALAXY IDENTIFICATION AND 
STATISTICS 

The galaxy finding procedure used here was a modified den- 
sity peak prescription: Based on the results of KKK91, KK, 
cooled particles are preferred tracers of the galaxy distribu- 
tion. We therefore began by searching for local maxima of 
the density field of cold particles exceeding a predetermined 
threshold value Sth- (The local maxima were also required 
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Figure 8. Slices of BSI-25* simulation at z=0. a) unshocked; b) 
cooled; c) 10 5 K < T < 1.5 X 10 6 K; d) T > 10 7 K. 



to be maxima within a 5 cell neighborhood.) The mass of 
the corresponding galaxy was then computed by summing 
the contributions of all particles within 0.5 cell lengths of 
the maximum in the lower resolution simulations and within 
one cell lengths of the maximum in the high resolution sim- 
ulations. This leads to the same physical length for selecting 
halos in simulations of the same box size, i.e. to comparable 
catalogues. Since the centroid can shift due to addition of 
new particles, the procedure was iterated six times, and a 
centroid, velocity, and mass were assigned to the galaxy. For 
each simulation, several galaxy catalogs were constructed us- 
ing various choices of threshold. Appropriate choices lie in 
the range S t h ~ (1.5— 3)as- In the end, one seeks a catalog in- 
cluding at least as many galaxies as would be expected from 
observational counts of galaxies with luminosities exceed- 
ing the characteristic Schechter luminosity L, in the given 
volume. For numerical simulations, the mass corresponding 
to "L»" -galaxies can only be determined from this require- 
ment a posteriori. We observed that reducing the threshold 
simply resulted in including additional galaxies at the low 
end of the mass spectrum, but had virtually no effect on 
the galaxies at higher masses. Thus we can extract one big 
galaxy catalog from each simulation, which can be cut off at 
the low-mass end as required. A natural cut-off is implied 
by the finite resolution of the simulations. At later stages, 
this cut-off grows due to the overmerging effect. 

In Fig. H we compare mass functions derived from the 
small box simulations where we expect the most reliable 
identification of single galactic halos. The halo masses span 
the range from 10 10 M Q to 4 x 10 12 M Q , and they can be 
fitted by Schechter curves of the form 

dn/dlogM = n»(M/M,)" p exp(-M/M»), (6) 
For BSI-25* we get M* = 10 12 M Q , n» = 2.6 x 10~ 2 /i 3 Mpc" 3 
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Figure 9. The 'galaxy' mass function of BSI-25 (thin solid his- 
togram), of BSI-25* (thick solid histogram) compared with the 
CDM-25 function (dashed histogram). The solid and dash-dotted 
curves correspond to comparison Schechter curves with parame- 
ters provided in the text. 
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Figure 11. The cumulative halo mass function of BSI-200* at 
redshifts z = 1.8,0.67,0.43,0.25,0.1,0. The mass range corre- 
sponds to small groups of galaxies which undergo substantial 
growth after z = 1.8. 
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Figure 10. The halo mass function of BSI-25 at redshifts z = 
2, 1.5, 1, 0.5, 0. For z = a best fit Schechter distribution is shown 
using a characteristic mass M t = 6 X 10 11 Mq, a power index 
p = —1.1 and mass density n* = 1.4 X 10 -2 /i 3 Mpc -3 . 

and p — 0.8, compare the solid line in Fig. ^. The lower 
dashed line represents the observational results from Efs- 
tathiou etal. (1988), n, = 1.6 x 10~ 2 /i 3 Mpc -3 , p = 0.1, 
and also M* = 1O 12 M0 (assuming a quite high M/L ratio 
of about IOOM0/L0). Similar parameters of the Schechter 
function are derived in Loveday etal. (1992). 

The simulated mass spectra lead to approximately the 
same abundance of M* galaxies as the observations. Com- 
parison of galaxy mass spectra between CDM and BSI sim- 
ulations at the same resolution shows that the CDM leads 
to a slower inclination at small masses, despite of the ex- 
cessive power at small scales. Due to overmerging (Katz & 
White 1993; Kauffmann etal. 1993), however, we know that 
spectra from PM simulations may not be directly compared 
with observational data. One way to imagine how overmerg- 
ing comes about is as follows: Potential wells associated with 
galaxies can generally bind particles with velocities of about 
250 km s -1 . However, the potential wells in clusters can bind 



particles with velocities of about 1000 km s _1 . The proba- 
bility that a particle within a cluster will be near one of the 
deeper troughs may be slightly elevated, but not enough to 
distinguish galaxies, at least not using the present scheme at 
the presently attainable resolution. In a hydrodynamic sim- 
ulation with proper treatment of cooling (Katz & Weinberg, 
1992, Katz & White, 1993, and Evrard, Summers & Davis 
1994), one would expect baryons to condense earlier in the 
galaxy (halo) potential wells, before the galaxies are assem- 
bled in the cluster. On the other hand, the PM dynamics 
leads to an merging of halos per se. We do not intend to 
make predictions on the galaxy mass functions in BSI from 
our simulations. We only use the halo identification scheme 
to get a catalog of objects with a reasonable number density, 
whose clustering properties should be characteristic for the 
primordial perturbation spectrum. 

In order to understand the development of (over-) merg- 
ing in a qualitative way, we studied the evolution of the mass 
spectra of halos (massive galaxies and groups) . The redshift 
dependence of the differential galaxy mass function of BSI- 
25 is shown in Fig. [l(j. A Schechter fit is possible for all 
redshifts, only at z — the number of very massive halos 
begins to exceed the exponential fall-off. Note that while 
the number of low-mass halos (M < 1O U M ) grows by less 
than a factor of two between z = 1 and z = 0, the number of 
galaxies of M > 5 x 1O 11 M0 grows by a factor of about 10 or 
more. The redshift dependence of cumulative mass functions 
for BSI-200* is shown in Fig. [TlJ Focusing our attention on 
the behavior of the curves beyond 2 x 1O 13 M , we found 
an increasing of the number of high-mass halos after about 
z — 0.67, suggesting that at earlier epochs (before cluster 
formation), overmerging may be less severe. 

Summarizing, CDM appears to lead to more massive 
objects than BSI (compare in Fig. |^ the thin solid and 
dashed lines for BSI and CDM distributions, respectively). 
However, in view of overmerging, the apparent disagreement 
of mass spectra derived in both scenarios with mass spectra 
inferred from luminosity functions cannot be interpreted as 
a weakness of either model at the present stage. 

We can gain important information from the spatial 
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Figure 12. The correlation function of galactic halos (identified 
in BSI-200*, (dash-dotted line) and of cluster halos (using BSI- 
500, solid line). For comparison we show a straight line with slope 
1.6 

distribution and velocities of galaxies that are derived from 
our catalogs. 

In Fig. [l^, we show the two-point correlation function 
for the galaxy catalog constructed from the simulation BSI- 
200*. It is well described by a power law £ = (r/ro)~ 7 for 
lft _1 Mpc < r < 15/i^" 1 Mpc, with slope 7 w 1.6 and corre- 
lation length ro ~ 6ft~ Mpc. Both values are in reasonable 
agreement with galaxy surveys. (Because the existence of a 
power law for the slope of the correlation function is a quite 
stable property of hierarchical clustering, as demonstrated 
in numerous simulations, it is not a strong discriminator be- 
tween models.) The slope is a bit on the shallow side, indi- 
cating that our spectrum could have a slight power deficit at 
galaxy scales. The correlation radius is also bit high; the ha- 
los in this galaxy catalog may be biased toward high masses 
and probably include small groups. Nevertheless, consider- 
ing the uncertainties in identifying galaxies and the effect 
of overmerging, the BSI simulation results are quite com- 
patible with observations. On the contrary, the correlation 
function for 'galaxies' identified in the COBE-normalized 
CDM simulations is much too steep (and one might expect 
that without overmerging it would have been even steeper). 
Finally, the break in the correlation function at 20ft _1 Mpc 
and the zero crossing over 25ft _1 Mpc are quite remarkable. 

A robust measure of clustering is provided by the count- 
in-cell variances, which are volume averages of the two-point 
correlation function. The (mass- weighted) variance of galaxy 
counts in spheres of radius 8h~ Mpc in the BSI-200 simula- 
tion is consistent with the observed variance 1. We compare 
the variances of galaxy and dark matter counts for spheres 
of different radii to get a measure of the integral bias, 

6 M = ^galM/^DtvlM- (7) 

To minimize the effect of overmerging, we use in any case 
mass-weighted cell variances. Fig. 13 shows this bias param- 
eter as a function of the chosen length scale. For the CDM 
models with more power on galactic scales, we obtain a slight 
antibiasing, while the bias of the BSI halos varies from 1 to 
2.5, if one inspects the simulations BSI-25 to BSI-500. It 
should be noted that the bias factors within one simulation 



Figure 13. Integral bias of different BSI vs. CDM simulations, 
the error bar are 1 — cr errors estimated from the count in cell 
variances. 

are almost constant. For deriving the counts-in-cells, we used 
randomly placed spheres with a total volume not exceeding 
the sample volume. Therefore we obtain statistically signif- 
icant results. The mean errors in the cell counts are used to 
estimate the error bars in the 'bias' in Fig. 13. The some- 
what surprising antibias of the 'peak-selected' galaxies in 
CDM is connected with the low threshold of the galaxy cat- 
alog used in this analysis (compare the mass functions of 
the galaxies) . On the other hand this result is in accordance 
with the normalization of the corresponding simulations, the 
comparison of the linear theory a dm (8ft _1 Mpc) variances, 
cp. Fig. 4, and the measured unit variance of galaxy counts. 
The results of the integral bias are summarized in Table 2. 



8 ANALYSIS OF GALAXY CLUSTERS 

Clusters of galaxies are widely acknowledged as a sensitive 
test of cosmological scenarios. Their formation is connected 
with the very recent decoupling of large masses from cos- 
mic evolution. Henry and Arnaud (1991) were the first to 
derive restrictions on the power index of the primordial per- 
turbation spectrum from X-ray cluster data. Using a cat- 
alog of 25 clusters with X-ray flux Lx > 3 x 10 -11 erg 
cm _2 s _1 (flux measured in the range 2 — 10 keV), they de- 
rived the dependence of the spatial number density on the 
X-ray temperature. Since the latter can be used as an in- 
dicator of of the cluster mass (cp. the estimates of Evrard 
1990), they obtained an estimate of the condensation prob- 
ability of different mass scales. The results of Henry and Ar- 
naud (1991) are based on the Press-Schechter (1974) theory, 
which supposes a Gaussian distribution of rare mass fluctu- 
ations when smoothed with a top-hat filter of the respective 
scale. Henry & Arnaud used in their analysis a power law 
spectrum P(k) oc k n from which they obtained a spectral 
index n ~ -1.7 (-2.4 < n < -1.3) at k ~ 0.3ft Mpc" 1 . This 
is shallower than the standard CDM model on the relevant 
scale (compare also Bartlett & Silk, 1993). Here we identify 
galaxy clusters in the largest simulations (BSI/CDM-500) 
using a effective search radius of 1.5ft Mpc, the Abell ra- 
dius. The resulting mass function of cluster halos is shown in 
Fig. 0. The observational points stem from a collection of 
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Table 2. Integral bias of galaxy halos 



Simulation 


radii of spheres r 
h _1 Mpc 


mean bias 
{b) 


standard 
deviation 


BSI-200* 


5 


- 40 


1. 1 


0.2 


BSI-500 


12.5 


- 100 


2.5 


0.5 


BSI-200 


5 


- 40 


1.7 


0.1 


BSI-75 


2 


- 12.5 


1.0 


0.1 


BSI-25 


0.8 


5 


O.S 


0.3 


CDM-500 


12.5 


- 100 


0.9 


0.3 


CDM-200 


5 


- 40 


0.8 


0.2 


CDM-75 


2 


- 12.5 


0.6 


0.1 


CDM-25 


O.S 


5 


0.5 


0.2 




M [M ] 



Figure 14. The cumulative number density of cluster halos of 
the BSI model (solid line, BSI-500) is a good description of the 
data of Cen & Bahcall (1993); full dots: optical data mainly from 
Abell clusters; squares: X-ray observations from EXOSAT and 
HEAO II. CDM simulations (dashed line) significantly overpro- 
duce cluster halos. 

both optical and X-ray data superimposed to get an univer- 
sal empirical mass function (Cen & Bahcall 1993). The BSI- 
simulation fits the observational data quite well, although 
at M < 6 x 1O 13 M we see a slight overproduction of clus- 
ters. Cluster formation is a quite recent process, as Fig. |ll| 
demonstrates. The CDM model leads to a tremendous over- 
production of halos; in particular, the mass function is too 
steep on the cluster scale, this conclusion is typical for Q — 1 
models. The break length of BSI is strongly restricted by the 
cluster mass function as shown by a Press-Schechter analy- 
sis (Miiller 1994a). 

The cluster number density is a sensitive test of the 
power spectrum at the transition scale between the differ- 
ent inflationary stages. An independent test of the power 
spectrum on these scales is the cluster-cluster correlation 
function (Bahcall & Soneira 1983, Klypin & Kopylov 1983). 
While the BSI correlation function, shown in Fig. |l2| as a 
dashed line, has the correct slope, it has a comparatively 
small correlation radius, ro ~ 10 ft _1 Mpc. Only slightly 
larger correlation radii are derived by Dalton et al. using the 
APM survey (r = 14.3 ±2.35 /i _1 Mpc) and by Collins etal. 
(1994) for a ROSAT selected cluster survey. Our estimate of 
ro ~ 10/i -1 Mpc may be marginal consistent with the opti- 
cal catalog. Of similar importance is the range of the region 



of positive correlation of cluster-cluster correlation. We find 
from our simulation a positive cluster-cluster correlation up 
to at least 80/i~ 1 Mpc. The same seems to be indicated by 
(up till now somewhat uncertain) observational results. 



9 LARGE AND SMALL-SCALE VELOCITY 
FIELDS 

Using catalogs of objects (large galaxies or clusters) con- 
structed from the 200/i _1 Mpc (M > 2 x 10 12 A'/ Q ) and 
500/i _1 Mpc (M > 3 x 10 13 M Q ) simulations, we estimated 
typical bulk velocities (and standard deviations) by com- 
puting unweighted average velocities in spheres of radius R 
(from 2.5 to 75/i _1 Mpc) centered about 512 equally spaced 
observers. A top-hat filter was used. (Note that with 8 3 
spheres considered there is of course some overlap in the 
larger spheres.) The results of the 500 /i _1 Mpc simulation 
are shown in Fig. [if]. Bertschinger etal. (1990) estimated 
the average velocities of the mass within spheres of radius 
40 and 60/i _1 Mpc centered at the Local Group as 388 ± 67 
and 327±82 km s _1 , respectively (see data points in Fig. |l4]) . 
At the 1-cr level, these estimates are compatible both with 
the BSI results and with CDM. The rms velocities for BSI- 
200* (not plotted) and BSI-500 are both well fitted by curves 
of the form w rms = vo exp(— x/xq), with vo ~ 490 km s , 
x = 41.7/i _1 Mpc for simulation BSI-200*, and v fa 520 
km s _1 , xq = 66.7/t _1 Mpc for simulation BSI-500. The min- 
imum mass included in the catalogs was 2.6 x 10 n M Q for the 
high-resolution simulation BSI-200*, 2 x 10 11 M for the sim- 
ulations BSI-200 and CDM-200, and 3 x 10 13 A'/ Q for the sim- 
ulations BSI-500 and CDM-500. However, varying the mini- 
mum mass had little effect (i.e., resulted in a shift well within 
the range of statistical uncertainty) . The Lauer and Postman 
(1994) measurement (689 ± 178) km s _1 at 150/i _1 Mpc lies 
above the bulk velocity in any of 512 bins for all CDM and 
BSI simulations at the largest radius computed (75/i _1 Mpc). 

An important property of BSI models as studied in 
GMS94 is to provide a natural mechanism for a relatively 
high Mach-number or "cold" flow (Ostriker & Suto 1990, 
Strauss et al., 1993). Small-scale velocity fields may be stud- 
ied by measuring the rms line-of-sight relative peculiar ve- 
locity derived from galaxy catalogs (Turner 1976; Davis & 
Peebles 1983; Mo etal. 1993). 

Davis and Peebles (1983) first inferred the small-scale 
velocity field of galaxy clustering from the anisotropy of the 
small-scale velocity dispersion. They derived the rms relative 
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r [/iT'Mpc] 

Figure 15. Streaming motions for BSI-500 (lower solid line) and 
CDM— 500 (lower dashed line) for the largest simulation box. The 
thin solid line and dashed line delimit the upper one-sigma range 
obtained from different observers in the BSI-500 and CDM-500 
simulation, respectively. The two error bars apply to the data 
from Bcrtschingcr etal. 1990). 



radial peculiar velocities «rms — (i v i —V2) 2 ) of galaxy pairs 
(where v± and t>2 are the radial peculiar velocities) as a func- 
tion of the projected distance r v between them on the sky. 
In the simulations we derive the relative peculiar streaming 
motion as function of the projected distance in bins up to 4 
h~ Mpc, and the rms of the relative peculiar velocity differ- 
ences for galaxies with projected separation between 1 and 
2 ft -1 Mpc, compare Table 3. While the streaming motion 
falls off with increasing r, the velocity dispersion urms is al- 
most constant, therefore we give only one value. To restrict 
the contamination due to foreground and background galax- 
ies, we consider only pairs of galaxies with real separation 
less than 10 h~ 1 Mpc. The values of the Table 3 are directly 
comparable to the Davis and Peebles (1983) results which 
give rms peculiar velocities from the apparent anisotropy 
of the spatial correlation function using the CfA redshift 
survey: v r ms(R) = (340 ± 4Q)(r p h _1 Mpc) - 13 km s" 1 , for 
r < l/i -1 Mpc. We get somewhat smaller velocity differences 
for all BSI simulations, which result from the combined ef- 
fect of restricted power in the simulation box, and the insuf- 
ficient resolution of small mass elements. In particular, we 
cannot resolve single galaxies in the large simulation boxes, 
but only infer the peculiar velocities of 'groups' or 'clus- 
ters'. The highest peculiar velocities are predicted for the 75 
/i _1 Mpc and the high resolution 400 /i -1 Mpc simulations. 
On the contrary, the velocity dispersion of CDM galaxies is 
too high. 

Recent studies of the same and newer data sample seem 
to indicate that the small-scale velocity dispersion depends 
on the selection of galaxy types, especially whether large 
clusters are included or not (Mo etal., 1993). In the latter 
case, a significantly higher velocity dispersion is obtained. 
Also studying the groups and clusters in simulations, one 
gets higher rms velocities. These results do not refer to the 
analysis of the mean of relative velocities of all close pro- 
jected pairs as done here. 



10 CONCLUSIONS 

A substantial body of independent and plausible observa- 
tional evidence indicates that the primordial perturbation 
spectrum is not well described by the standard CDM model 
with Harrison-Zel'dovich primordial fluctuation spectrum. 
Here, we have seen that the two-parameter spectrum with 
broken scale invariance resulting from a double-inflationary 
scenario provides a good description to a substantial range 
of observational data: the galaxy number density, galaxy- 
galaxy correlations, the count-in-cell statistics of galaxies, 
the cluster mass and X-ray temperature functions, cluster- 
cluster correlations, small and large-scale velocity fields. The 
observed filamentary structure is evident in the galaxy dis- 
tribution, the distribution of cold particles and (most promi- 
nently) in the distribution of warm (~ 10 6 K) gas. 

With regard to counts in cells, the very encouraging 
result — obtained from the distribution of DM particles - 
is that the BSI simulations fit the data quite well in a wide 
range of scales, [(5 — 60)/i _1 Mpc]. In contrast, the CDM 
counts in cells are inconsistent with observations at about 
the 2-cr level. 

The variance of galaxy counts depends strongly on the 
chosen galaxy identification scheme. The results are sum- 
marized in Table 2 with the estimated integral bias. For 
simulations BSI-200 and BSI-200*, the results are at best 
typical for galactic halos, while in BSI-500 we measure in- 
stead cluster halos. In contrast, for the CDM simulations we 
found 'galaxies' that are slightly anti-biased with respect to 
the dark-matter distribution. 

BSI seems to produce fewer massive objects than CDM. 
However, in view of overmerging, care must be taken in inter- 
preting the apparent disagreement of mass spectra derived 
in either scenario with mass spectra inferred from luminosity 
functions. 

The slope and correlation length of the BSI galaxy cor- 
relation function as discussed in Section 7 are close to ob- 
served values. The cluster correlation function has a corre- 
lation radius of (10 — 15)/i~ Mpc and stays positive at least 
out to about 80/i _1 Mpc, consistent with our observational 
knowledge. 

Rms line-of-sight velocities from BSI simulations appear 
to be closer to the Davis & Peebles (1983) data than compa- 
rable CDM simulations, which lie systematically too high as 
mentioned already in the first papers on CDM simulations. 

Large-scale streaming velocities in the 500/i _1 Mpc sim- 
ulations are well within 1-cr statistical variations as deter- 
mined from the simulations. Considering the current uncer- 
tainties in these measurements, they pose no problems for 
the model at this time. Moreover, a circa 25 per cent in- 
crease in the normalization as suggested by the Gorski et al. 
(1994) reanalysis of the COBE data would suggest roughly 
a 7 per cent boost in typical large-scale motions at the scales 
in question. But in any case our predictions are much lower 
than the Lauer & Postman (1994) results, our BSI expecta- 
tion at a scale of 150 /i -1 Mpc is about 55 km/s. 

An important test for the BSI spectrum will be the more 
precise measurement of multipole moments of the microwave 
anisotropy at scales of a few hundred Mpc (the Doppler 
peak). The advantage of these measurements is that they 
directly sample the linear part of the spectrum. A decisive 
question for BSI will be a more precise observational deter- 
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Table 3. First (f s tr; defined positive inwards) and second (urms) moments of distribution 
of line-of-sight projected relative peculiar velocity differences of galaxy halo pairs. 



Simulation 


v s fa [km s 1 ] 
(0-1) 
h~ 1 Mpc 


(1 - 2) 
7i _1 Mpc 


(2 - 3) 
/i _1 Mpc 


(3 - 4) 
/i _1 Mpc 


vims [km s x ] 
(1-2) 
h~ 1 Mpc 


RQT 900* 


183 


159 


142 


130 


256 


BSI-25* 


50 


40 


27 


18 


192 


BSI-500 


176 


169 


139 


106 


243 


BSI-200 


141 


117 


111 


101 


241 


BSI-75 


106 


104 


97 


88 


253 


BSI-25 


82 


78 


69 


62 


192 


CDM-500 


362 


373 


316 


266 


469 


CDM-200 


224 


216 


201 


172 


484 


CDM-75 


204 


218 


207 


199 


644 


CDM-25 


100 


126 


131 


110 


469 



mination of the epoch of formation of the first objects. As 
mentioned in Sect. 5, improved treatment of ionization and 
heating processes could allow use of quasar abundance and 
quasar absorption line data to obtain more stringent limits 
on structure evolution in BSI. Of comparable importance 
are data expected to be available in the near future on the 
evolution of the cluster mass function at redshifts z « 0.5. 
This problem presents a theoretical challenge, because it re- 
quires an accurate hydrodynamical treatment of cluster gas 
dynamics in the context of large-scale structure simulations. 

The most important success of the BSI scenario is the 
excellent agreement with present observational data on very 
large scale structure. In this paper, we have seen through 
numerical simulations that smaller scales are generally in 
excellent agreement with available data. In the future, it 
will be important to investigate additional characteristics of 
large-scale structure using methods such as topological stud- 
ies, percolation analysis and computation of void probability 
functions. 

In order to study mass spectra and statistics in a more 
reliable manner, hydrodynamic simulations should be per- 
formed. However, the results of our PM simulations strongly 
suggest that the BSI model studied here deserves full at- 
tention as a reliable description of cosmological structure 
formation. 
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